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Abstract 

Extensions of GMSB models were recently explored to explain the recent reports 
of the Higgs boson mass around 124 — 126 GeV. Some models predict a large \i 
term, which can spoil the vacuum stability of the universe. We study two GMSB 
extensions: i) the model with a large trilinear coupling of the top squark, and ii) 
that with extra vector-like matters. In both models, the vacuum stability condition 
provides upper bounds on the gluino mass if combined with the muon g — 2. The 
whole parameter region is expected to be covered by LHC at y/s = 14TeV. The 
analysis is also applied to the mSUGRA models with the vector-like matters. 



1 Introduction 



Gauge mediated SUSY breaking (GMSB) models are one of the best motivated extensions 
of the Standard Model (SM) at the TeV scale. For instance, dangerous flavor-changing 
neutral currents are naturally suppressed, and the big-bang nucleosynthesis can work 
successfully without suffering from photo- and hadrodissociation by the gravitino decay. 

Recently, the ATLAS and CMS collaborations have reported excesses of events in the 
SM Higgs boson searches at the mass around 124 — 126 GeV [I]. In the minimal SUSY 
Standard Model (MSSM), such a relatively heavy Higgs boson mass requires either a large 
trilinear coupling of the top squark or very large squark masses. In GMSB models, since 
the trilinear coupling of the top squark is predicted to be small, a large SUSY-breaking 
mass scale is required to explain such a heavy Higgs boson. 

On the other hand, the experimental result of the muon anomalous magnetic moment 
(muon g — 2) [2 J is deviated from the SM prediction at more than the 3a level [3]. The 
deviation may be explained by SUSY contributions, which are required to satisfy 

(V(SUSY) = (26.1 ± 8.0) x 1(T 10 (1) 

at the la level. However, if the soft mass scale is large, SUSY contributions to the muon 
g — 2 become suppressed. As a result, it is very difficult to explain both the Higgs boson 
mass of 124 — 126 GeV and the muon g — 2 anomaly in GMSB models. 

Recently, extensions of the GMSB models were studied in light of the Higgs boson 
mass and the muon g — 2 [U EJ [6J [7] . Some of them predict a large Higgsino mass term 
(/i term). Interestingly, large [i helps to enhance the SUSY contributions to the muon 
g — 2, since the smuon-Bino diagram is proportional to the fi parameter [8]. However, 
such a large \i parameter can destabilize the vacuum of the universe, since the trilinear 
coupling of the stau-Higgs term is enhanced (see [9j [10J E] for recent studies). In this 
work, the vacuum stability bound is studied in the following two GMSB extensions: i) 
the model with a large trilinear coupling of the top squark [121 E], and ii) that with extra 
vector-like matters [U 131 [5] . In both models, the trilinear coupling of the stau becomes 
large in the region where the Higgs is heavy. It will be found that the stability condition 
provides upper bounds on the gluino mass if combined with the muon g — 2. We will also 
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show that the whole parameter region which is favored by the Higgs mass and the muon 
— 2 is expected to be covered by the LHC experiment at y/s = 14 TeV. The vacuum 
stability bound on the mSUGRA models with the vector-like matters [H [H] will also be 
discussed, paying attention to the dark matter abundance. 

2 Vacuum Stability Bound 

The stability of the "ordinary" vacuum, i.e., the electroweak-breaking vacuum in our 
universe, may be spoiled, if the trilinear coupling of the stau, 

C ~ l^utan/3^/1 + h.c, (2) 

is too large. Here we omit the sub-leading contributions which are not enhanced by the 
Higgsino mass parameter, fi, nor a ratio of the Higgs VEVs, tan 0. As the trilinear 
coupling increases, an electric charge-breaking minimum becomes deeper. The stau 
trilinear coupling is therefore bounded from above so that the lifetime of the ordinary 
vacuum is longer than the age of the universe. 

The transition rate of the metastable vacuum is estimated by a semiclassical method, 
searching for so-called bounce solutions [15]. In Ref. [11], an approximate formula for the 
bound on /i tan /3 is obtained by using multi-dimensional bounce configurations, including 
top-stop radiative corrections to the Higgs potential: 

//tan/3 < 76.9y/M^m~ + 38.7(m fL +m? R ) - 1.04 x 10 4 GeV, (3) 

where m? L and rrif R are soft scalar mass parameters for the left- and right-handed staus, 
respectively. Although the correction to the Higgs potential depends on the masses of the 
stops, the bound (j3J) is affected only at the percent level [11] 

The above result is obtained in the limit of the zero temperature. The vacuum can 
transit through thermal effects in the early universe. The thermal decay rate of the false 
vacuum is usually estimated by following the method in Ref. [TB]. Evaluating the Higgs 
potential at the one-loop level, which includes the thermal potential coming from the top 



x The analysis in Ref. [11] assumed the A-term in the trilinear coupling of the top squark, At, vanishing. 
Although a large A t can contribute through the top-stop loops, the following conclusion is considered to 
be less sensitive to it, because the corrections from the top-stop loops are small in Ref. [TTj . 
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quark and electroweak gauge bosons, the stability bound on ([2]) can become more severe 
than the zero-temperature result ([3]) for a small stau mass region [10J. However, the 
bound on the stau trilinear coupling fl2]) would be more severe only up to ~ 10% [170], and 
hence we will adopt the bound in the following models, which provides a conservative 
constraint on the parameter space. 

The ordinary vacuum may be required to be absolutely stable if the vacuum 
expectation values of the scalar fields stayed away from the ordinary one in the early 
universe. In the following analysis, we will also show the constraint on the parameter 
space which is obtained by requiring the absolute stability of the vacuum (bearing in 
mind that the bound depends on the cosmological history). The upper bound on 
becomes more severe by about 50% compared to the metastable one ([3]) [IT] . 

Finally, let us touch on another possibility of the vacuum instability. In a class of 
SUSY models, the trilinear coupling of the top squark is predicted to be large in order to 
raise the Higgs mass. Such a large trilinear coupling may spoil the vacuum stability in the 
stop-Higgs plane similarly to the stau-Higgs plane discussed above. The metastability 
bound has been studied in Ref. [U] and obtained as A\ + 3/i 2 < 7.5(m^ +m~ R ), where A t 
is a SUSY-braking trilinear coupling of the top squark defined as V = —y t A t H®t* L t R + h.c.. 
On the other hand, depending on the thermal history of the universe, the vacuum may 
be required to be absolute stable. The condition is found in Ref. [IB] , which we will use 
in the following analysis. 

3 GMSB with Large A t 

In the simplest GMSB models, since the trilinear coupling of the top squark is suppressed 
at the messenger scale, it is difficult to realize the Higgs mass of > 124 GeV in the region 
where the muon g — 2 discrepancy is explained. This difficulty can be ameliorated if 
the Higgs doublets mix with the messenger doublets [121 E] ■ This mixing induces a large 
trilinear coupling of the top squark, and hence the Higgs mass is enhanced. In this section, 
we study the vacuum stability bound in this model. 

The model has six input parameters, (A mess , M mess , iV mess , tan/5, y' t , sgn(/z)), where 
A mess = F/M mess determines the soft mass scale. The sign of fi is fixed to be positive 
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in the following analysis, since otherwise the muon g — 2 anomaly cannot be explained. 
The Yukawa coupling y' t parametrizes the new superpotential, 

W = y^ L Q L 3f R , (4) 

where $£ is the leptonic component of the messenger field, and Ql% and Tr are the 
left- and right-handed top quark chiral multiplets, respectively. This coupling induces 
corrections for the soft parameters. The detailed mass spectrum is obtained in Refs. [T2|[7]. 
Remarkably, the soft SUSY breaking parameters of the third generation squarks, including 
At and m~ , receive extra threshold corrections. Furthermore, the soft mass of the up- 
type Higgs gets a negative contribution, which essentially enhances \i. 

The Higgs mass and the muon g — 2 in this model have been studied in Ref. [7J. It was 
found that the two conditions, > 124 GeV and a A1 (SUSY) > 10~ 9 , are simultaneously 
satisfied when the messenger scale is low, i.e., x ~ 0.4 — 0.5 with x = A mess /M mess , and 
the extra Yukawa coupling is large such as y[ > 1 at the messenger scale. It was also 
argued that large tan /3 is favored to explain the muon g — 2 anomaly. Since \x becomes 
quite large, the vacuum stability bound in this model is of great importance as discussed 
in the previous section. 

In Fig. [U contours of the Higgs mass and the muon g — 2 are displayed as a function of 
the physical gluino mass and tan/3. The other parameters are fixed to be (N mess ,x,y' t ) = 
(1,0.4,1.3), (2,0.45,1.2) and (3,0.5,1.4), respectively in each plot. For iV mcss = 2 and 
3, it is assumed that a pair of the messenger fields has a Yukawa interaction with the 
top. The formulae of the soft parameters can be found in Refs. [12], [7]. In the numerical 
analysis, the mass spectrum is obtained by modifying SoftSUSY [TH], and the Higgs mass 
and the muon g — 2 are evaluated by FeynHiggs [20]. In the plots, the green region denotes 
124 GeV < rrih < 126 GeV, and the muon — 2 is explained in the orange (yellow) region 
within the l(2)a level. It is found that the green band covers a wide parameter region 
because large y' t contributes to A t . On the other hand, the soft mass scale, i.e., the gluino 
mass, is constrained by the muon g — 2, and the region depends on tan/3. 

The metastability bound is shown by the black solid lines in Fig. [TJ This provides a 
severe upper bound on tan f3. Combined with the muon g — 2, the gluino mass is bounded 
from above by the stability condition. This result is rather stable against variations of 
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the model parameters, x and y' t . Around the parameter region where the Higgs mass is 
about 124 GeV, the SUSY contributions to the muon g — 2 do not depend much on x and 
y' t , though the Higgs mass is sensitive to them (see Ref. [7] for the dependence on x and 
y' t ). As a result, the gluino mass is found to be bounded as m g < 1.3 (0.9), 1.8 (1.4) and 
2.1 (1.5) TeV for N mess = 1,2 and 3, respectively, in the region where the muon (7 — 2 is 
explained at 2a (la) level. 

The absolute stable bound is also displayed by the black dashed line as labelled by 
f in Fig. [TJ If the thermal history of the universe requires the absolute stability of the 
ordinary vacuum, the constraint on tan/3 is found to be very stringent. As a result, the 
gluino mass becomes tightly limited such as as m g < 1.1 (0.8), 1.5 (1.2) and 1.7 (1.3) TeV 
for iV mess = 1, 2 and 3, respectively, for the muon g — 2 explanation at 2a (la). The 
stability constraint gives severe bounds on the superparticle masses, which is essential for 
the collider searches for the superparticles. We will discuss LHC searches later. 

Let us mention some other constraints. In the figures, the trilinear coupling of the 
top squark is large due to sizable y' t , which may destabilize the vacuum in the stop- 
Higgs plane. We have checked that the metastability condition is satisfied in the whole 
parameter region in Fig. [TJ On the other hand, the absolute stability bound is shown by 
the black dashed line with a label, t, for N mess = 2 in Fig. [TJ Above the line, there is a 
deeper vacuum than the ordinary one. Thus, most of the parameter region is potentially 
excluded by the constraint in this case, though whether it applies or not depends on the 
thermal history. 

A large trilinear coupling of the top squark also enhances a SUSY contribution to the 
branching ratio of b — > sj. The experimental result, Br(B — > X s 7) exp = (3.55 ± 0.24 ± 
0.09) x 10~ 4 [21], agrees with the SM prediction, Br(B ->> X s7 ) SM = (3.15 ± 0.23) x 
10 -4 |22j . very well, where B represents B° or B~ . Hence, the SUSY contributions 
are limited. We evaluated the SUSY contributions of the model at the NLO level by 
SusyBSG [23J. In addition to the experimental and SM uncertainties, errors of 10% are 
taken into account both for the SUSY and charged Higgs contributions, respectively (see 
e.g. [23] )• In the minimal flavor violation, the dominant contribution of SUSY comes from 
the chargino-stop loop diagram and is enhanced by large tan /3 as well as A t . Requiring 
that the SUSY contribution is within the 2a range of the SM and experimental results, 
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b — > S7 provides a constraint, which is displayed by the gray shaded region in the panel 
of iV mess = 2 and that with small tan (5 for iV m ess — 3 in Fig. [TJ 

3.1 LHC phenomenology 

The LHC constraints and prospects of the model have been discussed in Ref. [7]. Here, we 
discuss the LHC phenomenology in light of the vacuum stability bound obtained above. 
We will also extend the analysis including the cases of prompt NLSP decays. 

The collider signature of the SUSY events is determined by the next-to-lightest 
superparticle (NLSP), since the lightest superparticle (LSP) is the gravitino. Above the 
light blue dashed line the stau is the NLSP in the panels of N mess = 1 and 2 in Fig. (TJ while 
the neutralino becomes lighter below it. In the case of N mess = 3, the stau is the NLSP 
in the whole parameter region. The gray dot-dashed contours denote the stau mass. As 
tan (3 increases, the lightest stau mass decreases. The gray shaded region in the iV mess = 1 
panel and that in a large tan /3 region for N mess = 3 are excluded by a tachyonic stau. 

The collider signature depends on the lifetime of the NLSP. Its decay length depends 
on the mass of the gravitino and NLSP as cr ~ 10/im(m 3 / 2 /l eV) 2 (mNLSp/100 GeV)~ 5 . 
The gravitino mass is determined by the magnitude of the total SUSY breaking scale, 
Ftotai, which is equal to or larger than the F of the messenger sector, F = A mess M mess . 
As discussed in Ref. [7], the vacuum stability for the messenger sector tends to require 
F <C F tota i, and thus a long-lived NLSP is favored. For completeness, we also discuss 
the case of F ~ i^totai> where the NLSP can decay promptly q In the whole parameter 
region of Fig. m F/(V3M P ) ~ (1 - 10) eV, where M P is the reduced Plank mass scale. 
The NLSP either decays promptly or is long-lived, depending on the gravitino mass. In 
the following we discuss SUSY searches in the LHC for both cases. In the following, the 
Monte Carlo simulation relies on Pythia [21] to obtain the event shape, and the detector 



simulation is based on PGS [22]. The cross sections are evaluated by Prospino [2S] at 
the NLO and LO levels for the colored superparticle productions and the electroweak 
processes, respectively. 

When the messenger number is N mess = 1, the Bino is the NLSP in the whole parameter 



2 If the gravitino mass is mildly light, the NLSP can decay inside the detectors. Then, SUSY events 
include signatures with kinks or non-pointing photons. Such cases need to be studied elsewhere. 
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region where the metastability bound is satisfied, i.e., below the black solid line in Fig. [TJ 
Since the Bino is a neutral particle, the SUSY events are associated by a large missing 
transverse energy, E^ 1SS , if the Bino is long-lived. Thus, the hadronic channels with large 
missing energy are suitable for the detection, and the relevant production channels are 
those of the gluino and squarks, pp —> gg, gq and qq. Applying the high mass cut used 
in the ATLAS search [27], it is found that the gluino mass is required to be larger than 
around 900 GeV in Fig. [TJ Thus a corner of the parameter space is left to explain the 
muon g — 2 anomaly at the 1 a level. 

It is commented that this search is similar to those for the CMSSM models J27J [281 
I2"9"l I3T)] . The difference is that the sleptons are likely to be much lighter in the GMSB 
models. If the slepton mass is close to that of the NLSP, the leptons which are generated 
in decay chains of the colored superparticles are likely to be soft. Thus, the number of 
events which pass through the cut condition, the lepton veto, increases compared to that 
of CMSSM. 

If the Bino decays promptly, photons are likely to be generated when the Bino 
decays into the gravitino. Since the gravitino provides a missing energy, the di-photon 
with large E™ 1SS is a promising channel for the selection. Such a signature has been 
studied in ATLAS [3 lj and CMS [32J. The main production channels are electroweak 
productions of the charginos, neutralinos and sleptons. Applying the same cut as the 
ATLAS analysis [31], the bound is obtained as m g > 1.1 TeV for iV mess = 1 in Fig. [TJ 
Thus, if the Bino is the NLSP with a short lifetime, the LHC has already excluded the 
region which is allowed by the vacuum stability and is consistent with the muon g — 2 at 
the la level. 

As the messenger number increases, the stau tends to be lighter than the Bino. For 
-^mess = 2 and 3, the stau becomes the NLSP in most of the relevant parameter region as 
shown in Fig. [TJ If the gravitino is much heavier than 0(1) eV, the stau is long-lived and 
leaves a charged track in the detectors. Heavy stable charged particles are searched for 
in ATLAS [33] and CMS [31] . The CMS analysis was updated very recently by using the 
integrated luminosity, 4.7fb _1 [35]. If the updated analysis is done as the previous one [33] , 
the constraint is found to be rrtg > 1.9 TeV both for A^ mess = 2 and 3 in Fig. [TJ where the 
main production channels are electroweak productions of the charginos, neutralinos and 
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sleptons. The bound on the stau mass by the stau direct production is weaker than the 
ino productions. As a result, most of the parameter space which is favored by the muon 
g — 2 and the vacuum stability has been already excluded by the long-lived stau searches. 

If the stau decays promptly, the most relevant channel is associated by the di-lepton 
with large i?™ lss . The bound has been studied in the framework of the general gauge 
mediation, and the same-sign di-lepton channel provides the severest bound (36]. The 
ATLAS analyzed the opposite-sign di-lepton channel in the GMSB models [37J. In 
both cases, the gluino mass > 1 TeV is large enough to avoid the detections. Thus, 
no constraints are obtained for N mess = 2 and 3 in Fig. [1] as long as the stau is the NLSP 
and decays promptly. 

Finally, let us mention future prospects for the discovery. When the NLSP Bino 
is long-lived, the signal is very similar to the CMSSM. Since the gluino and squark 
masses determine the production cross section, the discovery potential can be read off 
by translating these masses into (m ,mi/2)- If the LHC will run at yfs = 7 TeV during 
2012, the gluino mass region of > 1.2 TeV is expected to be accessed for the integrated 
luminosity of 30 fb -1 [38]. The sensitivity would be better if the energy is upgraded to be 
y/s = 8 TeV. All the region favored by the muon (7 — 2 is possibly covered by the LHC at 
y/s = 14 TeV with the integrated luminosity of 0(1 — 10) ftT 1 [391 HQ]. On the other hand, 
if the stau is the NLSP and decays promptly, the multi-lepton channels are considered 
to be suitable discovery channels. Based on the study of the discovery potential of the 
low-scale GMSB, the gluino with a mass of around 2 TeV is expected to be discovered for 
0(1 - 10) fb" 1 at y/s = 14 TeV [H]. Thus, the whole parameter region which is favored 
by the muon g — 2 and constrained by the metastability bound can be covered by the 
future LHC experiments. 

4 The Model with Vector-like Matters 

We consider an extended model of the MSSM with vector-like matters. The extra matters 
are introduced as complete SU(5) multiplets, 10 = (O/, U', E') and 10 = (Q', U', E'). The 
details are found, e.g. in [4]. The model is highly motivated since it can explain the muon 
g — 2 and the Higgs mass, rrih — 125 GeV, simultaneously both for the GMSB [3J [5] and 
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mSUGRA scenarios [HE]. The superpotential is given by 

W = W MSSM + Y'Q'HJJ' + Y"Q'H d U' + M Q ,Q'Q' + M^U'U' + M E ,E'E\ (5) 

where Wmssm is the superpotential for the MSSM matters. In the presence of the extra 
up- type (s) quarks, the Higgs mass is lifted by their radiative corrections. This becomes 
significant when the Yukawa coupling is as large as 1" ~ 1, and the SUSY invariant 
masses, Mqi and M[//, are in the TeV scale [321 03j 03]- A model was proposed in 
Ref. [45J, which can naturally explain the existence of such TeV scale vector-like matters. 
The extra Yukawa coupling, Y', has an infrared fix point at Y'(EW) ~ 1 (33]- 

The vector-like matters raise the fi parameter. When there are extra matters, the 
soft mass squared of the up-type Higgs, rn 2 Hu , is driven down at the electroweak scale 
compared to that of the MSSM. This is due to the renormalization group (RG) evolution 
of m 2 Hu . The additional contribution from the extra squarks to the RG equation is given 
by 

(16vr 2 )^| 9 6Y' 2 (mk + m Q' + <> + A*) (6) 

at one-loop level, where Q is a renormalization scale, itiqi and niyi are soft masses of 
the extra squarks, and A' is a scalar trilinear coupling which corresponds to the second 
term of the right-hand side of Eq. (jSJ). Since Y' ~ 1 is required to lift up the Higgs 
mass, this correction becomes sizable. From the electroweak symmetry breaking (EWSB) 
condition, /i is found to be large, which can induce a deep electric charge-breaking 
minimum depending on tan /3 and the stau masses as discussed in Sec. [2j 

It is commented that the trilinear coupling of the top squark is small in the extra 
matter models. This is because large Y' suppresses A t during the RG evolution. Thus, 
A t at the weak scale is insensitive to the input at a high scale and has an infrared fixed 
point at a small value (331 E]- This feature is distinct from the model in Sec. El where 
large A t potentially spoils the absolute stability condition and can suffer from the b — > 57 
bound. The models with the vector-like matters are free from these constraints. 

In the following two subsections, the Higgs mass and the muon g—2 are evaluated under 
the stability condition both in the GMSB and mSUGRA scenarios. The mass spectrum 
of the SUSY particles are evaluated by SuSpect [46J, which is modified to include the 
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contributions to renormalization group equations from the vector-like matters at the two- 
loop level. The MSSM contributions to the Higgs mass is calculated by FeynHiggs at 
the NLO level, and those from the vector-like matters are at the one-loop level. The 
muon g — 2 is also calculated by FeynHiggs. In the mSUGRA, we calculate the relic 
abundance of the lightest neutralino by micrOMEGAs [47J. The collider Monte Carlo is 
simulated by the same packages as those in Sec. [3J Pythia, PGS and Prospino. In the 
following analysis, we assume Y" ~ 0, since otherwise the extra down-type squark reduces 
the SUSY contributions to the Higgs mass (see Ref. [5]). 

4.1 GMSB models with vector-like matters 

First we consider the simplest GMSB model, which is parametrized by A mess , M mess and 
N mess as well as tan/3 and sgn(/i). Here we set the messenger number as N mess = 1, since 
larger iV mess spoils the perturbative gauge coupling unification, unless the messenger scale 
is high enough. 

In Fig. |2l the Higgs mass, the muon g — 2 and the stability bounds are shown in the 
m^-tan plane for various messenger scales. The green band shows the Higgs mass of 
124 GeV < m° h < 126 GeV for Mq> = M v > = 0.6 and 1.0 TeV. In the orange (yellow) 
region, the muon g — 2 discrepancy is explained at the la (2a) level. 

The vacuum stability gives stringent constraints on the soft mass scale. Above the 
black solid line, the region is excluded by the metastability condition 03]). If the bound 
is combined with the muon g — 2, the soft mass scale, which is represented by the gluino 
mass in Fig. 121 is constrained from above. From ([3]), the upper bound is obtained as 
nig < 1.2 (1.7) TeV for the muon g — 2 at the la (2a) level when the messenger scale 
is 10 6 GeV. Also, in the region above the black dot-dashed line, the ordinary vacuum is 
not absolutely stable but becomes metastable. This gives stronger upper bound on tan f3. 
The bound is found to be m g < 1.0 (1.3) TeV for the muon g — 2 at the la (2a) level and 
M mess = 10 6 GeV. For both constraints, the upper bound on the gluino mass becomes 
severer if M mess is larger, 

In Fig. [21 the NLSP is the lightest stau (neutralino) above (below) the blue dashed 
line. In all the panels of Fig. [21 the NLSP is likely to be long-lived, since M mess is 
larger than 10 6 GeV. When the stau is the NLSP, the parameter regions below the black 
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solid line and above the blue dashed one seem to survive the metastability condition 
especially for low M mess . However, those regions have been already excluded by the 
searches for the heavy stable charged particle in the LHC as studied in Ref. [5]. On the 
other hand, if the messenger scale is as low as ~ 10 5 GeV, the gravitino mass could be of 
0(1 — 10) eV, and thus, the stau decays promptly. Then direct search bounds from the 
LHC are weak, and thus, the region with nig > 1 TeV is still available (see [361 137]). Since 
the gluino mass is tightly bounded by the metastability condition, the whole regions are 
expected to be covered by searches for the multi-lepton events in early stage of the LHC 
at y/s = 14 TeV [U]. The regions below the black dot-dashed line already excluded the 
stau NLSP. 

When the NLSP is the neutralino, long-lived neutralinos contribute to the missing 
transverse energy at the colliders. If the high mass cut of the ATLAS search [27] is 
applied, it is found that the gluino mass is required to be larger than around 750 GeV in 
all the panels of Fig. [2j Thus, a large part of the parameter space is still valid against 
the collider searches. According to the discussion in [5] , the whole parameter region with 
m g < 1.7 TeV is expected to be covered in future at the LHC if the integrated luminosity 
reaches O(10 - 100) fb" 1 at ^/s = 14 TeV. 

If the gravitino mass is 0(1 — 10) eV, the neutralino NLSP decays promptly in the 
detectors. The di-photon events with a large missing transverse energy are searched for 
and provide a constraint on the SUSY particle masses. If the ATLAS result [31] is applied, 
the bound is obtained as m g > 1.2 TeV for M mess = 2 x 10 5 GeV, in which case the mass 
spectrum is almost same as those of M mess = 10 6 GeV in Fig. |5J Thus, most of the 
parameter region which is favored by the muon g — 2 at the la level has been excluded, 
and the whole la region of the muon (7 — 2 is expected to be accessible in the LHC soon 
at v/i= 14 TeV [H]. 

4.2 mSUGRA models with vector-like matters 

Next, we consider the mSUGRA boundary condition, which is characterized by the 
following five parameters: the universal soft scalar mass tuq, the universal gaugino mass 
mi/2, the universal trilinear coupling A , tan/3 and sgn(/i) at the GUT scale. Since A t 
is suppressed at the weak scale and the other A^s are irrelevant for the Higgs mass and 
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the muon g — 2, we set A = in numerical analysis. We also take sgn(/x) = +1, which 
is favored by the muon g — 2. 

The contours of the Higgs mass and the muon g — 2 are shown in Fig. [3] and Fig. HJ 
where the definitions of the regions and lines are same as Fig. In Fig. [31 tan /3 are varied 
among the panels, while mo is set to be zero in Fig. HI The metastability bound from the 
stau is represented by the black solid line. This constraint becomes severer as tan /3 is 
larger. In addition, the region with larger m-o avoids the constraint. The latter feature is 
understood as follows: the stau mass is approximately proportional to m , while the size 
of fj, is less sensitive to it and rather determined by the gaugino (gluino) mass m.1/2 (see 
[H]) in most of the parameter region of our interest. Thus, the bound (j5J) is satisfied in 
a larger m region. The absolute stability bound displayed by the black dot-dashed line 
is understood similarly. 

From Fig. [31 it is found that when tan /3 is large and mo is small, the muon g — 2 can 
be of order 10~ 9 easily in the region where the Higgs mass is 124 — 126 GeV. However, 
such regions are likely to be excluded by the metastability condition. Moreover, if the 
ordinary vacuum is required to be absolutely stable, the lower bound on m as well as 
the upper bound on tan (3 becomes tighter. As a result, it is more difficult to satisfy both 
rrih = 124 — 126 GeV and explain the muon g — 2 anomaly at the 2a level unless tan /3 is 
properly small. 

In Fig. HI we show the Higgs mass and the muon g — 2 on the m^-tan (3 plane. Here 
mo is set to be zero so that the SUSY contributions to the muon g — 2 become as large 
as possible. Note also that the physical gluino mass is smaller by O(10)% than mi/2 in 
the extra vector-like matter models |44j . It is found that the gluino mass is bounded 
from above if the stability bounds on tan/3 are combined with the muon g — 2. The 
metastability bound is found to be m~ g < 0.8 (1.1) TeV when the muon g — 2 deviation is 
explained at the la (2a) level. If the absolute stability is imposed, the gluino is required 
to satisfy rrig < 600 (850) GeV for the la (2a) level of the muon g — 2. 

In the mSUGRA models, the Bino is the LSP. The searches for the multi-jets events 
with a large missing transverse energy provide a constraint for the gluino and squark 
masses. In the current model, the squark mass is related to the gluino one, and the gluino 
mass is required to be larger than around 750 GeV if the high mass cut of the ATLAS 
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search [27] is applied. Thus, most of the region which is favored by the muon g — 2 at the 
la level has been excluded already, and is expected to be covered by the LHC soon [38] . 
Moreover, since the vector-like matter is required to be relatively light in order to raise 
the Higgs mass sufficiently, the extra matters are expected to be discovered by searches 
for the fourth generation quarks at the LHC [5]. On the other hand, if the muon g — 2 
allows the 2a discrepancy, the upper bound on the gluino mass is relaxed. The gluino 
mass of nig = l.lTeV is in reach of the LHC experiment at i/i = 7 — 8TeV with the 
integrated luminosity of ~ 10fb _1 [38J. Thus, the whole favored region with the muon 
g — 2 at the 2a level will be checked in the near future at the LHC. 

Finally let us mention the relic abundance of the lightest neutralino. The blue solid 
lines denote a contour of the neutralino abundance with Ocdm^ 2 — 0.1 in Fig. [3] and Fig. HJ 
It is noticed that the region exists for small mo if tan (3 is large, where the coannihilation 
works efficiently for the neutralino LSP. Such parameters are likely to spoil the stability 
of the vacuum as mentioned above. Consequently, the coannihilation region is found for 
rrig < 600 GeV from Fig. HI However, the region has been already excluded by the LHC. 
Thus, in order for the neutralino LSP to explain the measured abundance of the dark 
matter, late-time entropy productions should occur after the neutralino decoupled from 
the thermal bath. 

5 Conclusion 

The simple GMSB models are disfavored by the muon g — 2 in light of the recent results of 
the Higgs boson search. Extensions of the models were recently explored to accommodate 
the muon g — 2 to the Higgs mass of 124 — 126 GeV. As a characteristic feature, some 
models predict a large \x term. In this work, the vacuum stability constraint from the stau 
was studied. In particular, we studied two GMSB extensions: i) the model with a large 
trilinear coupling of the top squark, and ii) that with the vector-like matters. In both 
models, the vacuum stability condition was found to provide upper bounds on the gluino 
mass if the constraint is combined with the muon g — 2. Based on the gluino mass upper 
bound, the discoveries at the LHC were also discussed. It was concluded that the whole 
parameter region obtained by the metastability bound and the muon g — 2 is expected to 
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be covered by LHC at y/s = 14 TeV. 

This analysis was also applied to the mSUGRA models with the vector-like matters. 
The whole region was found to be accessible in the near future at the LHC. It was also 
argued that the coannihilation region has been already excluded by the LHC experiments. 
Thus, dilutions are required to take place in early universe, since otherwise the relic 
neutralino overdoses the universe, or the R-parity might be violated slightly. 
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Figure 1: Contours of the Higgs boson mass of 124 — 126 GeV (green region) and the 
muon g — 2 at the l(2)er level (orange (yellow) region) in the GMSB model with a large 
trilinear coupling of the top squark. The sign of \i is set to be positive. The black solid 
and dashed lines labelled by f represent the upper bounds on tan/3 obtained from the 
meta- and absolute stabilities of the vacuum along the stau-Higgs direction, respectively. 
The region above the black solid lines are excluded. For A^ mess = 2, the black dashed 
line labelled by t corresponds to the upper bound on tan (5 obtained from the absolute 
stability of the vacuum along the stop-Higgs direction. For iV mess = 1 and 2, above the 
light blue dashed line, the stau is the NLSP, while the neutralino is lighter below it. For 
-^mess = 2 and 3, the stau mass is displayed by gray dot-dashed lines. The gray shaded 
regions are excluded by the tachyonic stau and b —> 57. 
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Figure 2: Contours of the Higgs mass, the muon g — 2 for various messenger scales in the 
GMSB model with extra vector-like matters. The region consistent with 124 GeV < rrih < 
126 GeV is shown as the green band for Mqi{= Mu>) = 600 GeV and 1 TeV. The orange 
(yellow) region is consistent with the muon g — 2 at the la (2a) level. On the blue dashed 
line, masses of the lightest neutralino and stau are equal to each other. Above (bellow) 
the line, the lightest stau (neutralino) is the NLSP. Above the black dot dashed line, the 
ordinary vacuum is metastable. The region above the black solid line is excluded by the 
vacuum instability. In the gray shaded region the stau mass becomes < 100 GeV. 
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Figure 3: Contours of the Higgs mass, the muon g — 2 and the relic abundance of the dark 
matter in mSUGRA models with extra vector-like matters. The region/lines are same as 
in Fig. [2] except for the blue solid line, which denotes the contour of the relic abundance 
of the lightest neutralino, f2cDM^ 2 — 0.1. In the left region of the black dot dashed line, 
the ordinary vacuum is metastable. In the left region of the black solid line is excluded 
by the vacuum stability condition. The other mSUGRA parameters are taken as Aq = 
and sgn(/i)=l. The gray region is excluded, since the stau is tachyonic or the chargino is 
too light. 
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Figure 4: The regions/lines are same as in Fig. [3j Here, mi/2 and tan/3 are varied with 
setting m = 0, A = and sign(//) = l. The region above the black solid line is excluded. 
The gray region is excluded, since the stau is tachyonic. 
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